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Abstract 

In this paper we calculate the gravitational-wave power-spectra 
corresponding to different types of transitions between the inflationary 
regime and the radiation era. We study four cases, where inflation is 
followed by a stiff matter phase, or by a dust dominated phase, or by 
a combination of the two, before the universe enters the radiation era. 
Use is made of differential equations for the Bogoliubov coefficients, 
allowing us to model all the transitions between the different eras as 
continuous. New features appear, for frequencies above 10 7 rad s _1 in 
an otherwise flat part of the spectrum. 

I. Introduction 

Although gravitational-waves of cosmological origin have not yet been 
directly detected, they are at present the object of an important effort of 
research, as they provide us with a unique telescope to the very early stages 
of the formation of the universe. Using the upper bound set up by the COBE 
measurement of the total flux of energy coming from the metric perturba- 
tions P, we explore the consequences, to the power spectrum of the stochas- 
tic background of gravitational-waves, of the possible existence of different 
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regimes of transition between the inflationary universe and the standard ra- 
diation dominated expansion. This transition is far from being well known 
and it may be of some interest to explore the signatures associated with dif- 
ferent possibilities, just assuming that an ideal detector will be able, in the 
more or less distant future, to explore in some detail the power spectrum. 
Examples of this can be seen in the works of Giovannini ([2], [3] and [1]), as 
well as in Mendes and Liddle concerning models with thermal inflation [5]. 

In the present work, and contrary to these authors, we do not use the 
sudden transition approximation between different eras. We model the tran- 
sitions as continuous and integrate numerically the differential equations, 
derived by Parker [B], for the time dependent Bogoliubov coefficients. As- 
sociated with the sudden transition approximation we always have an over- 
production of gravitons of large frequencies, requiring the introduction of 
an explicit cut-off for frequencies larger than the rate of expansion at the 
time of the transition. Quite often, also, some more or less delicate cancel- 
lations have to be worked out in the limit of the low frequencies. All this is 
automatically taken into account by the use of continuous transitions. 

For simplicity, we assume that inflation is dominated by a cosmologi- 
cal constant, after which we have a transition period modelled by a self- 
interacting scalar field, whose energy is continuously transferred to a fluid, 
with an equation of state of the form p = jq, through the action of a fric- 
tional term. This is followed by the radiation and matter eras. At a certain 
point, during the matter era, the universe is supposed to become dominated 
by dark energy, described by an equation of state of the form pd e = WQd e , 
with w negative, in agreement with the results of the WMAP satellite]?]. 

In the next section we introduce our model and derive the main equations. 
In section III we describe the results of the numerical simulations for the four 
cases investigated and the corresponding spectra. In the first example, after 
inflation we have a period dominated by a stiff matter fluid; in the second 
example inflation is followed by a dust dominated period. In the other two 
examples, inflation is followed by a combination of stiff matter and dust. 
All the cases end with the universe entering the standard radiation era. In 
section IV we make some additional comments and discuss the possibility of 
detecting the effects described. 

II. The model 

The background metric is defined by 
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ds 2 = -a 2 (r])(-dr] 2 + 8 ij dx i dx i ) (1) 

At the beginning the expansion is driven by a cosmological constant, up 
to a time rj = rji, when both the transition period and our integration begin. 
The scale factor up to r\ = rji is given by 

where H is the Hubble constant during inflation and 771 = a(rji) j a! (rji) + rji 
[8]. In conformal time the Einstein's equation for the scale factor is 



where 



and 



a" = ^^a 3 (g - 3p) (3) 



Qtotai = ^{\f' 2 + a 2 V{^)) + q (4) 



Ptotai = ^(t^' 2 - a 2 V&)) + V (5) 

define the total energy density and pressure, and g and p = are the energy 
density and pressure of the fluid into which the scalar field energy density is 
transferred. For numerical convenience, we divide the integration into two 
parts, the first part taking us from rji up to rji, the end of the transition 
period, also marking the beginning of radiation domination. At this point, 
the energy that initially was in the scalar field will have been transferred to 
the radiation field p = g/3 (the energy contained in the matter and dark en- 
ergy components are still completely negligible). Once the present content of 
the universe is known from the WMAP measurements, we know the redshift 
corresponding to this point and, from it, we easily find the redshift corre- 
sponding to rj = rji. We then begin the second part of our integration that 
will take us from 771 till the present time, rj a , using the Einstein's equation in 
the form 

;/ AnG 3,g m o 1 i in Qdeo \ tn \ 

a = ( ^ + ( + 3|u;|) ^(iqH)) (6) 
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where, in conformal coordinates, the radiation fluid does not contribute to it 
and p mo and Qd eo are the present values of the energy densities in the matter 
and dark energy components. We define the scale factor at the present time 
to be a(r) ) = 1; knowing g mo and Qdeo, we know the redshifts corresponding 
to the transitions between radiation to matter domination and from matter 
to dark energy. 

The other background equations, for the first stage of integration, are the 
equation for the scalar field 

if" + 2%' + a 2 d v V{tp) = -al>' (7) 
a 

and the equation(s) for the fluid(s) p = jg. If, after the end of inflation at 
r) = rji, and before the radiation stage (7 = 1/3), we have an intermediate 
stage with p = jg (7 7^ 1/3), we then have the two equations 

g' = -3-{l + ^)g + -T v if' 2 -l.a.Tg (8) 

CL (X 

and 

g' = -4—Q r + aTg (9) 
a 

I being a control number, which is different from zero only after a pre-defined 
and more or less extended interval of time such that, by then, the energy of 
the scalar field has effectively been transferred to the fluid p = 7^; Y v and 
T are the frictional decay constants that control this tranfer of energy from 
the scalar field to the fluid and from this to radiation field, respectively. The 
remaining equations to be used in both stages of integration are the ones for 
the Bogoliubov coefficients, defined below (equations [[8] and [EH]) . We must 
stress that our aim is not to construct a truly realistic model for the transition 
period marking the end of inflation, but to explore different possibilities with 
the help of a simple and reasonable toy-model. 

As for the tensor perturbations (latin indices from 1 to 3) of the space 
part of the metric a{rf) 2 {5ij + /ij 3 (^,x)), they can be expanded, in the usual 
manner, in terms of plane-waves 



^[J (27r) 3 / 2 a(r/)V2k 



d 3 k 



x[a p (k, 77)£jj(k,p)e lk - x ^(k, rj) + herm.conj.} (10) 

where x denotes the spatial coordinates, k =| k |=27ra/A = cua; the index 
p runs over the two polarizations of the gravitational-waves, and is the 
polarization tensor, a p the annihilation operator and £ the mode function for 
the gravitational- waves. The mode function obeys the equation 

£" + (k 2 - a"/a)£ = 0, (11) 

the derivatives being with respect to conformal time rj. We now express the 
creation and annihilation operators, at rj, in terms of the initial creation and 
annihilation operators A],(k) and A p (k), through the Bogoliubov coefficients 
a(k, rj) and /3(k, rj): 

a(k, rj) = a(k, V )A(k) + (3*(k, r/)A f (k) (12) 
a and /3 satisfying the relation 

\a\ 2 - \(3\ 2 = 1. (13) 

These coefficients, which, in the sudden approximation, are calculated by 
requiring the mode functions and their derivatives to be continuous across 
the transitions, obey the set of coupled of differential equations [SI E] 

^(v) = ^Hr ] )+(3(r ] )e 2 ^-^ (14) 

and 

(3'( V ) = -^0%) + a( V )e- 2 ^^)- (15) 
2k a 

where r] is an arbitrary constant to be put equal to rji] in Parker's inte- 
gral equations the function W(k,rj) has here been replaced by the ansatz 
W(k,r]) = k. Introducing the functions X(k,r)) and Y(k, rf) through the 
definitions 

a = -{X + Y)e ik{ "-Vo) ( 16 ) 

13 = -(X -Y)e- ik ^~ 11o) (17) 
the two equations above take the form 
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a" 



X" + (k )X = (18) 

(X 

Y = % -X' (19) 
k 

these being the remaining equations of our system. When the scale factor 
obeys a simple power-law, the solution to equation (fTS!) is expressed in terms 
of Hankel functions. When this is not so, as in the present situation, we have 
to integrate it numerically, with appropriate initial conditions. For the case 
of an initial stage dominated by a cosmological constant and a scale factor a 
defined by equation (j2J), the solution is well-known and will define the initial 
conditions for X and X': 

X - 1 -H^) (20) 

= Y = 1 - - (21) 

k k 2 (rji - rji) 2 k(rj! - rji) 

with (rji — r)i) already known. At the end of the integration, when the initial 
and final states are characterized by different kinds of mode functions, an 
additional Bogoliubov transformation is necessary, as explained in references 
[S] and [H]; with a final stage dominated by dark energy with an equation 
of state pde = WQde, w ~ —0.78, we may ignore the final projection, as it is 
numerically irrelevant. 

Once we have X and X' at r\ = r] , a final and j3fi na i can be recovered from 
equations (TTdT) . (TTT1) and (1T91 . The number of gravitons created during the 
expansion of the universe is given by \(3f ina i\ 2 [6J; taking into account that 
the density of states is u 2 du/(27i 2 c 3 ) and that each graviton contributes with 
two polarization states 2hu>, then, from the definition of the energy density 
dE = P(u))dw, we have the following expression [10] for the power-spectrum 
P{u): 

P{") = ^Wnal? (22) 

in units erg s cm -3 . We can also express the results in terms of the relative 
logarithmic energy-spectrum of the gravitational-waves, which is defined by 
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1 dp gw 
pcdhiui 



(23) 



where p c is the value of the present time critical density and p gw is the 
gravitational-wave energy density 



When performing these calculations within the sudden transition approxi- 
mation, due to a logarithmic divergence, a cut-off needs to be introduced in 
equation (1241) on the high-frequency gravitons, which depends on the speed 
of the physical transition [10] . This is not necessary when the transitions are 
continuous, as shown in [8] and [11]. In the low-frequency limit we have a 
natural cut-off provided by the size of the present horizon [TD] . 

The gravitational-wave contribution cannot be larger than the total amount 
of anisotropy measured by COBE. This implies that Q is limited by (value 
taken from [2] ) 



for frequencies corresponding to the present horizon size a;^ or , where we used 
H = 71 km s _1 Mpc" 1 , giving Uh or = 10~ 17 rad s _1 . In the numerical 
simulations, we used an inflationary scale chosen in order to avoid any conflict 
with the known limits for the scalar perturbations. 

III. Numerical simulations 

To proceed with the numerical calculations, we used the simple scalar 
field potential 



taking p = 4 and \i = the initial value of (p. Thus, at the beginning of 

the integration V = V Q , V D being then equal to the scale of the inflationary 
regime. We now give and discuss a few examples. 

Example 1 

In the first example we assume that the inflationary regime is followed by 
a period dominated by a fluid with a stiff matter equation of state, equation 
P with 7 = 1. 




(24) 



Q < 1.37 x 10 



-10 



(25) 




(26) 
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An instance of a model of this kind is the quintessential inflationary model 
developed by Peebles and Vilenkin [T2] and explored, from the point of view 
of the production of gravitational waves, by Giovannini [3]. First, the energy 
in the scalar field is converted into a stiff matter fluid, by the action of 
the term; we consider the conversion concluded when QsuffCt 6 becomes 
constant. Then, after a period of time, which we took, in the present instance, 
approximately equal to Arj = 4.5 x 10 5 , in Planck units, the stiff matter is in 
turn converted into radiation with the help of the phenomenological frictional 
term V stiff in equations © and ©. We used as initial conditions, at r] — rji, 
the values ipt = 0.001, ip' { = —0.00001, = Tsuff = 0.01, all in Planck units; 
we also fixed the scale of inflation (Vq) 1 ^ = 4.3 x 10 16 GeV. Finally, from 
the results of the WMAP we fixed the equation of state of the dark energy 
component as p = —0.78g and fide = 0.73. 

With these values, the redshift corresponding to rji is Z{ n f = 8.53 x 10 28 . 
The relative logarithmic energy-spectrum Q(u) is given in figure 1. We see 
the new feature, first noticed by Giovannini [3] , appearing in the hard branch 
of Q(lj), for frequencies > 10 7 rad s -1 . From the results in [8], we expect a 
cut-off to emerge naturally for frequencies uj >z 10 12 rad s _1 , sensitive to the 
first transition between the end of inflation and the stiff matter stage. 

The jump in Q(u) not only is sensitive to the scale (K) 1 / 4 , as it also 
increases with the increase in the interval of time Arj, during which we allow 
the stiff matter stage to dominate. In the present case we have, for instance, 
n(uj hor ) = 1.3 x 10 -11 , = 10~ 12 rad s" 1 ) = 2.7 x 10" 13 , in the horizontal 
part of the spectrum, and Q(u = 10 11 rad s~ x ) = 9.6 x 10~ 13 , at the peak. 

For comparison we added the curve corresponding to the sudden transi- 
tion approximation, using equations (4.8) in [10J. Being a different model, 
we had to fix the redshift for the transition from the radiation to the matter 
era. We used the value z = 2.2 x 10 3 , midpoint between the values given in 
[7] for z eq and Zd ec - In each of the examples, figures 1-4, the Hubble constant 
at the sudden transition between inflation and radiation was chosen in or- 
der for Vt{ujhor) to have the same value as for the corresponding continuous 
transition case. This required a small adjustment in its value, explaining the 
slight difference in the horizontal line in the different plots. 

Example 2 

This time, following the inflationary period, we have a period with a 
dust-like equation of state p — 0. 

We keep the initial conditions of the preceeding example, except for = 
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0.01 and (K) 1/4 = 6.2 x 10 16 GeV, giving Q(uj hor ) = 5.6 x 10~ 13 . The 
extension of the matter dominated period (during which ^a 3 = constant) is 
Ar] = 1.4 x 10 5 . We have z- mi = 8. x 10 29 . 

Instead of an increase in Q, we find a decrease (see figure 2), which will also 
depend on Arj. In the present case £1(uj = 10 10 rad s _1 ) = 6.14 x 10~ 16 , lower 
than the value Q(u = 10~ 12 rad s" 1 ) = 6.7 x 10~ 15 . Although during the dust 
matter interval gravitons are produced, particularly during the transitions to 
and from this period, the increase in the scale factor a is faster than in the 
previous example, leading to a stronger redshift of the energy, and a decrease 
in fi. The increase around uo = 10 12 rad s _1 is due to the presence of the 
V((p) potential. 

In the same spirit of exploring different configurations between inflation 
and radiation, we next investigate what happens when we have a period with 
stiff matter followed by dust, or vice- versa, before the final transition into 
the radiation era. 

Example 3 

Inflation is now followed by a dust-like period (7 = 0), which is in turn 
followed by stiff matter (7 = 1), before the final transition into the radiation 
era. 

The dust and stiff matter periods have an extension given by Arjdust = 
1.8 x 10 5 and Ar] s tiff = 1.4 x 10 5 , respectively. The energy in <p is converted 
into dust with the help of T v , then dust is converted into stiff matter by 
the action of T d and, finally, stiff matter into radiation via F stiff] we keep 
T^ = T d = V a// = 0.01. The scale of inflation (V Q ) 1/4 = 7.6 x 10 16 GeV and 
Vt(uj hor ) = 7.2 x 10~ 13 ; z inf = 7.5 x 10 29 . 

The results are shown in figure 3. We notice that, in the range of fre- 
quencies uj > 10 7 rad s _1 , we first have an increase, due to the extra power 
brought in by the stiff matter component, and then the dip caused by the 
dust stage. It is illustrative to compare these results with those from example 
4. 

Example 4 

We exchange the two intermediate periods, first a stiff matter period and 
then dust. 

Taking Ar} stiff = 6. x 10 5 , Ar} d = 1.4 x 10 6 and (V Q ) 1/4 = 4.8 x 10 16 
GeV, with the same T's, we find the results seen in figure 4, first a dip and 
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then an increase with respect to the horizontal line, this last at the level of 
tt(uj = l(T 12 rad s" 1 ) = 4.1 x 1(T 13 . 

In the examples above we used values for (V^) 1 / 4 between 5 and 7.6 x 10 16 
GeV, in reasonable agreement with the approximate limit 2.6 x 10 16 GeV 
found in [13] . Given the uncertainties around this period, we believe these 
values for (K) 1 / 4 to be quite reasonable. 

Other, more complex, situations could have been devised, but these exam- 
ples already prove that, if we have a complicated transition between inflation 
and radiation, and an ideal detector, by observing the shape of Q,(u>) we 
may eventually be able to say a lot about such transitions, as the changes 
in Q(u>) not only reflect the different equations of state, but also the order 
and sequence of the different periods. This can be seen from the comparison 
between figures 3 and 4. In all cases studied, the size of the new features 
depends on the extension of the different intermediate periods and the scale 
taken for the inflationary period. 

IV. Concluding remarks 

In this paper we developed a simple toy model to study different types of 
transitions between the inflationary regime and the radiation era and, with 
the help of differential equations for the Bogoliubov coefficients, calculated 
the corresponding gravitational-wave power-spectra. The usual sudden tran- 
sition between inflation and radiation was replaced by a more or less complex 
period between these two regimes, period dominated by one or two fluids de- 
fined by equations of state of the form p = 7^ . All transitions involved were 
taken as continuous. Four examples were given, whose spectra can be seen 
in figures 1-4. 

Some additional comments for the stiff matter case, figure 1, are in order. 
A bump in the power-spectrum was found, beginning around 10 7 rad s _1 , 
but not a spike as big as the one reported in ref. 3. This is possibly due 
to our use of continuous transitions, automatically avoiding the inevitable 
over-production of high-frequency gravitons that take place in the sudden 
transition approximation. If we accept, for the scale of inflation, the limit 
derived in ref. 13 from WMAP, it is difficult to see how a much larger signal 
can be obtained. 

Concerning the possibilities of observing these new features, it is impor- 
tant to stress that it will not be enough to detect the flat part of the spectrum 
above 10 7 rad s -1 , we should also be able to detect the spectrum below this 
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frequency. Only then we will be certain that a new feature is indeed present. 
We shall have to wait for future generations of gravitational-wave detectors, 
in order to be able to see the main features described in this paper, which 
take place at frequencies above 10 7 rad s" 1 (see ref. 14, particularly page 
332). The same applies to the other examples we have studied. We agree 
with the final comment in ref. 3 that the GHz region should be carefully 
exploited. 
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Fig. 1 - The relative logarithmic energy-spectrum Q(u>) corresponding 
to example 1. We have inflation followed by a period dominated by a fluid 
described by a stiff matter equation of state. The thin curve shows the result 
for the sudden transition approximation, as defined in the text. In all figures 
we have the logarithms in the basis 10. 



12 



Log Q [to] 




Fig. 2 - The case described in example 2, where inflation is followed by a 
period dominated by a fluid with the dust equation of state p — 0. The thin 
line again refers to the sudden transition approximation. 
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Fig. 3 - Q for the case of example 3, where we have inflation followed by a 
period dominated by dust (p = 0), in turn followed by a fluid described by a 
stiff matter equation of state. The sudden transition results are represented 
by the thin line. 
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Fig. 4 - Q for the case of example 4, with inflation followed by stiff matter, 
in turn followed by dust. Thin line for the sudden transition approximation. 
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